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Abstract 

We report a discovery of steady long-living slowly eastward moving large-scale coherent twin 
cyclones, the equatorial modons, in the shallow water model in the equatorial beta-plane, the 
archetype model of the ocean and atmosphere dynamics in tropics. We start by constructing an¬ 
alytical asymptotic modon solutions in the non-divergent velocity approximation, and then show 
by simulations with a high-resolution numerical scheme that such configurations evolve into steady 
dipolar solutions of the full model. In the atmospheric context, the modons persist in the pres¬ 
ence of moist convection, being accompanied and enhanced by specific patterns of water-vapour 
condensation. 


* Corresponding author. Email: zeitlin@lmd.ens.fr 
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As is known, the large-scale atmosphere and ocean dynamics in tropics is largely de¬ 
termined by equatorial waves. Indeed, superposition of atmospheric Outgoing Longwave 
Radiation (OLR) data onto the spectrum of equatorial waves gives a very good correspon¬ 
dence [1], except for a branch corresponding to slow eastward-moving motions associated 
with the Madden-Julian Oscillation (MJO) [2], The Kelvin waves, the only species of the 
equatorial waves that can produce steady large-scale eastward-moving structures, have much 
faster velocity than MJO [3]. This fact motivates a search for coherent dynamical structures 
with slow eastward propagation which, thus, should be necessarily nonlinear. There do exist 
steady eastward-moving structures in a fluid layer on the mid-latitude /-plane, the modons 
which are exact dipolar solutions of the quasi-geostrophic (QG) equations [4], Although the 
classical modons were found analytically in the QG approximation, it was shown first by a 
computer-assisted analysis [5] and then by direct numerical simulations [6] that correspond¬ 
ing solutions exist also in the parent rotating shallow water (RSW) model on the /- plane. 
The QG modon solutions were extended to the full sphere [7, 8] in the framework of equiva¬ 
lent barotropic model, which is structurally close to QG. The pioneering paper [9] advanced 
an idea that the MJO could be related to such a modon in spherical geometry and showed 
that, like in the midlatitude tangent plane [6], the equivalent barotropic modon persists in 
the RSW on the sphere. Yet, the archetype model for understanding dynamics in the tropics 
is RSW in the equatorial beta-plane, e.g. [10]. We, thus, look for eastward-moving modons 
in this model. However, it encounters an obstacle from the very beginning, because there 
is no consistent QG approximation at the equator, e.g. [11], so the known modon solutions 
can not be borrowed. Our main observation is that there is a dynamical regime in RSW, 
called long-wave approximation in oceanography [12], which corresponds to small pressure 
variations, and gives, to the leading order, equations which do allow for modon solutions. 
In the atmospheric context this regime is the RSW analog of the non-divergent equatorial 
balance model [13]. The relevance of this model to large-scale tropical motions in the atmo¬ 
sphere is supported by scale and data analyses [14, 15]. We construct the modon solutions 
in this regime and, following [6], use these asymptotic solutions to initialize high-resolution 
numerical simulations with the full RSW model, and show that coherent dipolar steady 
eastward-moving structures do arise and persist. We then show that, in the atmospheric 
context, inclusion of moisture, with condensation and evaporation, enhances the modons, 
without disrupting long-time coherence, and produces specific convective patterns. 
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We start with RSW equations in the equatorial beta- plane with no dissipation: 


| dtv + v • Vv + f3 y z A v + gVh = 0, 

{ d t h + V • (vh) = 0 , 

where V = ( d x ,d y ), v = (u,v), u is zonal and v- meridional components of velocity, h is 
geopotential height (thickness), f3 is the meridional gradient of the Coriolis parameter, and z 
is a unit vertical vector. The interpretation of the model in the oceanic context is direct, while 
in the atmospheric context it should be understood as vertically averaged primitive equations 
in pseudo-height pressure coordinates [16], and can be extended to include water vapor with 
condensation and related latent heat release [17], as well as surface evaporation (see below). 
We introduce a pressure perturbation parameter A: h = H(1 + Xg), where H is unperturbed 
thickness, and fix spatial, velocity and time- scales: (x,y) ~ L, (u,v) ~ V, t ~ L/V. Under 
hypothesis that A —> 0, and gHX/V 2 = 0(1), hence V « y/gH, i.e. the characteristic 
velocity is much smaller than the phase velocity of the Kelvin waves c = \fgH, the non- 
dimensional equations take the form: 

d t v + v • Vv + /3y z A v + Vr? = 0, (2) 

X(d t ri + v • Vrj) + (1 + Xg) V • v = 0, (3) 

where /? = (3L 2 /V , and v = v 0 + Avi + .... In the leading order in A (3) gives V • v 0 = 0, 
the motion is non-divergent, and uq = —d v ip, vq = d x i[). Cross-differentiation of the zonal 
and meridional momentum equations results in the equation for the stream-function: 

V 2 r/>t + J ("0, V 2, 0) + f3xjp x = 0, (4) 


where J denotes the Jacobian. The modon solutions are built following [4], They are are 
obtained under hypothesis of steady motion with constant zonal velocity U , by supposing a 
linear relationship between the absolute vorticity and stream-function in co-moving frame, 
which gives inhomogeneous Helmholtz equation. It is solved by separation of variables in 
polar coordinates in terms of Bessel functions, first in the outer domain under condition of 
decay, and then in the inner domain, and matching the inner and outer solutions across a 
circle of a given radius a in the plane. The solution has the form: 


Vw = -K^g K i(P r ) sin #, r> a, 


'tpint 


Up 2 

k 2 J± (ka) 


Ji{kr) 


r 

k 2 


(! + [/ + Uk 2 ) sin0, r < a, 


( 5 ) 
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FIG. 1. Normalised streamlines and velocity field of an asymptotic modon in stationary (left ) and 
co-moving (right) frames. Dashed circle: separatrix of radius a. 


where J\ and K\ are ordinary and modified Bessel functions of order one, respectively, p 
is real, p 2 = /?/£/, so U > 0, and the motion is eastward. For each pair (a,p) there exist 
a series of eigenvalues k arising from matching conditions, of which the lowest corresponds 
to a dipole, to be called the asymptotic modon. The next eigenvalue gives a quadrupolar 
solution called “shielded modon”, and so on. By construction, this solution can be centred 
anywhere in the equatorial beta-plane, yet it is only when its center is at the equator that it 
results in a cyclonic pair. The streamlines of the asymptotic modon are presented in Fig. 1. 
First-order corrections tq, V\ can be found in the next order of the asymptotic expansion. 
There is no guarantee that such expansion converges, but numerical results presented below 
suggest that there is indeed an exact solution of (1) corresponding to this asymptotic one. 

We initialized the numerical simulations in the RSW model with the velocity field corre¬ 
sponding to the asymptotic modon solution (5) with U <C y/gH, and flat pressure. Simula¬ 
tions were performed with high-resolution well-balanced finite-volume numerical scheme [18] 
in a rectangular domain with sponges at the boundaries, which allow to mostly evacuate short 
inertia-gravity waves. No explicit dissipation was added. The domain was chosen to be sym¬ 
metric with respect to the equator and wide enough 9 L& x 6.5L^, where L& — (gH) 1 ^//3 1 / 2 is 
the equatorial deformation radius, in order to minimize the influence of boundaries onto the 
modons. We varied spatial resolution from 600 x 600 to 1200 x 1200, to check numerical con¬ 
vergence. The natural units of length and time in the numerical scheme are L^ and 1//3L^, 
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FIG. 2. Initial adjustment of the asymptotic modon with non-dimensional a = 0.5,17 = 0.15, as 
seen in the thickness and velocity fields. Time in units of l//3Ld- 


as respectively, which corresponds to nondimensional beta equal to one. Their numerical val- 

97 ues, with g and \i fixed for the Earth, are uniquely defined by the value of H. the equivalent 

98 depth. In the atmospheric context, for H ~ 10 km, Ld ~ 3000 km and 1 / f3Ld ~ 5 h. At the 

99 initial stages of the simulations, cf. Fig. 2 we observed an adjustment of initial configuration 

100 with emission of inertia-gravity waves and formation of a dipolar coherent structure, both 

101 in pressure and velocity fields. This dipolar structure in a form of twin cyclones then moves 
los eastward without changing form, as follows from Fig. 3. To exclude the influence of inertia- 
105 gravity waves produced by the initial adjustment and partially reflected by the boundaries, 
ioo the modon was “nudged”, i.e. the far gravity wave field was removed, and the simulation 


6 











































































- 4 - 3 - 2-101234 


X(L d ) 

FIG. 3. Snapshots of h for eastward-moving equatorial modon after initial adjustment at t = 
0,10, 20, 30,40, 50 [1 /(/3Ld)\, from top to bottom (nondimensional a = 0.5, U — 0 . 2 ). 


107 reinitialized without it. The detailed characteristics of the “exact” modon are presented in 
los Fig. 4. Its coherence can be inferred from the scatter plot of Bernoulli function vs potential 
loo vorticity (PV) in the upper-right panel of the Figure, which clearly gives a line. As is known, 
no e.g. [5], Bernoulli function and PV are functionally dependent for steady moving solutions of 
in RSW equations. The phase speed of the modon depends on initial a, U. The modons keep 
in moving eastward for the long time without losing their coherence. The modon’s velocity 
U 3 can be inferred from the distance it covered with respect to initial condition. The energy 
n 4 of the modon remains practically constant, while we detected a weak enstrophy loss, which 
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FIG. 4. Upper row: thickness h of the equatorial modon at t = 15 [1 / (left); scatter plot of 
Bernoulli function in the co-moving frame vs potential vorticity (PV) (right). Lower row. Zonal 
section of the equatorial modon across the center of each cyclone (left), and meridional section 
across the center of the modon (right), a = 0.5, U = 0.2. 


us diminishes with increasing resolution (several % loss at high resolution), not shown. The 
in difference between the asymptotic and “exact” modons is clear from the comparison of their 
ns relative vorticity in Figure 5, which shows that the “exact” modon is more compact and has 
HQ a larger peak vorticity. Although we did not explore in detail the space of parameters of the 
m modon solutions, we varied the parameter a, the modon’s radius. We observed that while 

122 for a < Ld equatorial modons always emerge, for a > L ( ], with the same velocity scale, the 

123 initial asymptotic modon transforms into a packet of equatorial Rossby waves moving west- 

124 ward. We also checked the robustness of the equatorial modon to the details of initialization, 

125 by either removing the outer field of the asymptotic modon, or misaligning its axis. In all 

8 






















126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 



-1 - 0.5 0 0.5 1 -1 - 0.5 0 0.5 1 


x (L d ) x T d ) 

FIG. 5. Relative vorticity of the asymptotic (left) vs “exact” (right) modons. “Exact” modon is 
more compact and more intense. 

such simulations the eastward-moving equatorial modon was always emerging after leaving 
some “debris” and/or wobbling at initial stages. This indicates that equatorial modons are 
attracting solutions. 

After having established existence of equatorial modons in the adiabatic environment, we 
switched on the effects of moisture within the diabatic atmospheric moist-convective RSW 

(mcRSW) [17]. The equations of the model read: 

/ 

dfV + v • Vv + /? r/ z A v = — g V h, 

< d t h + V-(yh) = -iP, (6) 

k d t Q + V-(Qv) = —P + E. 

Here Q is a bulk amount of water vapor in the air column, 7 is a parameter depending on 
the underlying stratification, P is the condensation sink, and E is the evaporation source of 
moisture, which are parameterized as follows: 

p= Q ~ Qa H(Q-Q a ), E = a\v\(Q s -Q)U(Q S -Q). (7) 

T 

Q s is a saturation moisture threshold, r is relaxation time, a is a parameter regulating 
evaporation, and %(...) denotes the Heaviside function. As in previous test simulations [17], 
we take a uniform initial moisture distribution Q t close to saturation: Qi = Q s — 0.01 with 
Q s = 0.9, a short relaxation time r equal to several time-steps of the numerical scheme, 
7 = 1 , and a = 0(1 0 -1 ). We performed exactly the same simulations as in the adiabatic 
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FIG. 6. Steady eastward propagation of the equatorial modon, as seen in PV field in moist- 
convective (left) and adiabatic (right) environments, with the same initial conditions, t = 
0,10, 20, 30,40 [l/(/3Ld)], from top to bottom, (a = 0.5, U = 0.2, a = 0.05). 


case. They showed that the equatorial modons keep their coherence in the moist-convective 
environment. Moreover, as is known, and confirmed in mcRSW [19, 20], condensation leads 
to intensification of the cyclonic vortices. Correspondingly, the potential vorticity (PV) 
anomalies become stronger, the size smaller, and the phase speed of the convectively-coupled 
modon becomes higher than that of its adiabatic counterpart. A comparison of evolution 
of the same modon in adiabatic and moist-convective environments is presented in Fig. 6. 
A typical pattern of the associated moisture field, and characteristic condensation patterns 
at the front and at the rear of the modon are shown in Fig. 7. Such patterns persist 
all along the simulation. We should stress that, although we initialized the simulations 
with non-divergent velocity field of the asymptotic modon, the “exact” modon develops and 
maintains a characteristic convergence/divergence pattern (not shown), which redistributes 
moisture in its core, and prompts condensation in specific zones. Notice that as condensation 
enhances cyclonic vorticity [19], the moist-convective modon intensifies, and moves faster 
than the dry one. 

To conclude, we established existence of long-living, slow eastward-moving dipolar coher¬ 
ent structures, the modons, in the shallow-water dynamics in the equatorial beta-plane. In 
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FIG. 7. (a = 0.5,7/ = 0.2, a = 0.05). Distribution of the water vapor (left) and superposition of 
pressure (black) and condensation (red) isopleths (right) at t = 30 [1 / (dL ( j j]. 


the atmospheric context, the modons maintain their coherence in the presence of moisture 
and are accompanied by a specific and robust moist convection pattern, if condensation 
and evaporation are present. The scales of these structures are conditioned by the inter¬ 
pretation of the RSW model. At given g and (3 the actual scales depend on the value of 
the unperturbed thickness H, or of the corresponding deformation radius L,i- Although we 
cannot make a firm link to MJO at this stage, some features of the equatorial modon do 
resemble those of dipolar structures associated with the MJO [21]. Work is in progress in 
this direction, in particular including improved parametrisation of moist processes [20] and 
vertical structure, which has specific features in MJO, cf. [22] and references therein. Inde¬ 
pendently of a possible link to MJO, the very existence of the eastward-moving equatorial 
modons, with their characteristic velocity pattern, slow eastward propagation, and relatively 
weak signature in pressure should be kept in mind in data analyses and simulations with 
“big” atmospheric and oceanic models. It is known [23] that long-wave pressure anoma¬ 
lies in RSW at the equator produce westward-propagating dipolar Rossby-wave packets and 
eastward-propagating Kelvin waves, in accordance with the Gill’s scenario [24]. The genesis 
of eastward-moving dipoles is an important topic, to be addressed elsewhere. Our analysis 
above suggests that comparable zonal and meridional scales of the initial disturbance, and 
relative smallness of pressure anomaly are necessary for that. 

Let us comment on the relation of equatorial modons we discovered to those known, or 
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hypothesized, in the literature. First, while the relation of the spherical modons to the 
standard QG ones on the midlatitude tangent plane was established in the limit of small 
spherical modon size [7], this proof, which uses the structural resemblance of equivalent 
barotropic equations on the sphere and QG equations, both incorporating a finite barotropic 
deformation radius, is not directly transposable to the equatorial beta-plane case. So the 
proof that our modons is a limiting case of the spherical ones is pending. Second, westward- 
propagating localized solutions, which were also called modons, were constructed in [25] as a 
generalization of westward-moving equatorial Rossby-wave solitons, and confirmed by direct 
numerical simulations in [23]. The long-wave scaling used in [25],[23] differs from ours, with 
the main difference residing in the disparity of zonal and meridional scales. A possibility of 
existence of eastward-propagating modons, if some ad hoc “heuristic” terms were added to 
the systematically derived asymptotic equations, was also evoked in [25]. Question remains 
whether our modons could be related to this conjecture, although the difference in scal¬ 
ings seems to prevent this. We should add, for completeness, that meridionally translating 
zonally-symmetric modon-like solutions in the vertical plane at the equator were derived in 
[26] in the context of the atmospheric equatorial boundary layer. 
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